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Abstract. We find the high energy asymptotics for the singular Weyl-Titchmarsh 
ra-functions and the associated spectral measures of perturbed spherical Schro- 
dinger operators (also known as Bessel operators). We apply this result to 
establish an improved local Borg-Marchenko theorem for Bessel operators as 
well as uniqueness theorems for the radial quantum scattering problem with 
nontrivial angular momentum. 



1. Introduction 

In this paper we will investigate perturbed spherical Schrodinger operators (also 
known as Bessel operators) 

d 2 1(1 + 1) . . , ^ 1 
ax x I 

where the potential q is real-valued satisfying 



(1.1) r = -— + V _ 9 ' +q{x), />--, x e K+ := (0,+oo), 



(12) a^T^ (R ) J X<1{X) Gil (0,l), l>~h 

[x(l-log(x))q(x) E L L (Q,1), l = -±. 

Note that we explicitly allow non-integer values of I such that we also cover the 
case of arbitrary space dimension n > 2, where 1(1 + 1) has to be replaced by 
1(1 + n — 2) + (n— l)(n — 3)/4 [38l Sec. 17. F], or the case of scattering of waves and 
particles in conical domains [Si. Due to its physical importance this equation has 
obtained much attention in the past and we refer for example to [T], [2], [TS], |28j . 
[52"] . [55] . [55] and the references therein. 

Of course one of the most interesting applications is the scattering problem in 
R 3 with a spherically symmetric potential. Originating in the seminal work of 
Heisenberg, the question if the potential q is determined by the scattering phase for 
one fixed value of the angular moment / was investigated by several authors. In the 
simplest case I = this question is completely understood by now and we refer, e.g, 
to the monographs [3], [TT], or [33]. However, for general I this question remained 
open to the best of our knowledge. As one of our main results we will show that, 
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in combination with the eigenvalues and norming constants, the scattering phase 
determines the potential uniquely. 

Our approach will be based on the Jost function [3D]. Basic properties of the 
Jost function in the case of general / were first formulated in |26j . Asymptotics 
for the Jost solutions and the Jost functions are well-known and were investigated 
by several authors [6], [16], [15]. Moreover, a scattering approach between the 
unperturbed operators for I — and I € (— ^, ^] was recently introduced in [2]. 

Our approach will use singular Weyl-Titchmarsh-Kodaira theory as one basic 
ingredient. This approach was first used by Kodaira [26] in connection with radial 
scattering theory as already mentioned above. This original work contained some 
gaps first pointed out by Kac |21j who proposed an alternate approach based on 
Krein's method of directing functionals. However, these results did not get much 
attention until recently when Gesztesy and Zinshenko |17j took it up again and 
triggered a large amount of results [5]. [5]. [12], [15 ] . (13 ) . [25 ] . [25 ] . [50 ] . [5T ] , 

Finally, the outline of our paper is as follows: We will first establish the high en- 
ergy asymptotics of the singular Weyl function and the spectral measure associated 
with operators of the type (|1.1[) . Then we will use this result in order to improve 
some inverse spectral uniqueness results for such operators from [8] and |28j . Fi- 
nally, we will apply our finding to the quantum mechanical scattering problem. 

2. Asymptotics of the Weyl function 

We will use r to describe the formal differential expression and H to describe 
the self-adjoint operator acting in i 2 (R + ) and given by r together with the usual 
boundary condition at x = 0: 

(2.1) ]imx l ((l + l)f(x)-xf(x)) = 0, fe [--A 

x— >0 Z A 

We are mainly interested in the case where t is limit point at oo, but if it is not, we 
simply choose another boundary condition there. Moreover, one could even replace 
M + by a bounded interval (0,6). 

Next, notice that our assumption (|1.2[) implies that the equation 

(2.2) ry = zy 

has an entire system of solutions (f>(z,x) and 9(z,x) such that 

(2.3, »(^)={§s g t*,x), !*:?:■ 

where </>, € W 1 ' 1 ^,!] and 0(0) = 0(0) = 1 (see, e.g., [31]). The singular Weyl 
function m : C \ M — > C is defined such that 

(2.4) ip(z,x) = 6(z,x) +m(z)(f>(z,x), A e C \ M, 

either belongs to L 2 (c, b) for some (and hence for all) c S (0, b) if H is limit point at 6, 
or satisfies the boundary condition at x = b if H is limit circle at b. Note that, while 
the first solution (j)(z,x) is unique under this normalization, the second solution 
0(z, x) is not, since for any entire f(z) the new solution 0(z, x) = 6(z, x) — f(z)<p(z, x) 
also satisfies (I2.3|) . Note that the corresponding singular m- function fh is given by 



(2.5) 



rh(z) = m(z) + f(z) 
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in this case. Moreover, it was shown in [28], [31] that the singular m- function (|2.4p 
admits the following integral representation 

(2.6) m(z) = g(z) + (1 + I ( -1^ - ^) Z * R. 

Here := [4 + f J ) the function <? is real entire, and p : R — > R is a nondecreasing 
function satisfying 

rt<rt p(A+)+p(A-) /• rfp(A) 

(2.7) p(A) = , p{0) = 0, / ,,,,,, < oo. 



(1 + A 2 )«' + ! 

The function p is called £/ie spectral function and e?p is me spectral measure. We 
also remark that the value of Ki is best possible (see [3T] for further details). 

In the special case I = the function m is the classical m-function and Marchenko 
[33) proved that in any nonreal sector in <C + the m-function satisfies 

(2.8) m(z) = — / = 2(l + o(l)), |z|^+oo, 

where the branch cut of the root is taken along the negative real axis. The estimate 
for the remainder term was later improved by Krein, Levitan, and Marchenko (see, 
e.g., |34) . [37] ). By now there is a vast literature on high energy asymptotics 
for m-functions of general Sturm-Liouville operators and it seems the first results 
were obtained by M.G. Krein and IS. Kac [22] [23l [24], Hille [19], Everitt [10] 
and Kasahara [25]. The most complete results on one-term asymptotics for the 
m-functions were obtained by Bennewitz [3] (cf. also our Appendix [A| . 

In the present paper we are interested in the asymptotic behavior of the singular 
m-function (|2.4p of perturbed Bessel operators. Our main result is the following 
extension of Marchenko's asymptotic formula (|2.8[l . 

Theorem 2.1. Suppose 

(2.9) H =-£ + l l£l + q{x) , 

where q satisfies (|1.2[) . Let m(-) be the singular m-function (|2.4[) . Then there is a 
real entire function g such that in any nonreal sector 

(2.10) m(z)-g(z)=mi{z)(l + o(l)), \z\ -> +oo, 
where 

zf 2 lHlo g (-z), i + I £ N„, r(! + |)2W 

Moreover, the spectral function satisfies 
(2.12) p(A)=ft(A)(l+o(l)), A^+c», 

where 



(2.11) m,(*) = ^ ^f+iW V ; 2 C, 



C* 2 ,.,.i±s . . 1 



(2-13) ft (A) = - 7 ^- 1 KX[o,oo)(A)A' + t, / > _1 

Our proof relies on two main ingredients: The one-term asymptotic formula for 
the m-function (Theorem IA.2P and the commutation methods for Bessel operators 
|30) . It will be given in Section [3] where we will use a Liouville type transform to 
establish a connection between the perturbed Bessel operators and the Krein string 
operators — r ^ d ^ ■ in Sectional we apply Theorem l2.1l to derive an improved local 
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Borg-Marchenko result for Bessel operators. In the final section [5j we apply our 
results to establish uniqueness for the quantum scattering problem. 

Necessary results on the asymptotic behavior of m-functions of regular Sturm- 
Liouville problems are collected in Appendix A. In Appendix [B] we provide some 
estimates and analytical properties of the solutions <p(z, x) and tj){z,x). 

3. Proof of Theorem 12.11 



First, let us note that it suffices to prove Theorem 12. II in the case when x = b is 
a regular endpoint. In fact, by |281 Lemma 7.1], we know 

(3.1) m( 2 ) = -5^4 + 0(- 1 



(Z,C) y^Z(t) 2 (z,c) 

as \z\ — > +oo in any nonreal sector. Hence the asymptotic behavior of m depends 
only on the behavior of q near and is independent of the behavior of q outside a 
neighborhood of as well as of a possible boundary condition at b. 

Therefore, without loss of generality we set 6 = 1 and assume that q G L x (c, 1) 
for some c G (0, 1). We divide the proof into three steps. 

3.1. Connection with the Krein string operator. In this section we restrict 
our considerations to the case when x = is in the limit circle case, that is 
I G [— |, 4). It is well-known that under this assumption the self-adjoint opera- 
tor H associated with r is lower semibounded. However, instead of the boundary 
condition (|2.1|) induced by <f> we will now use the one induced by 0{\q) plus an 
arbitrary boundary condition at x = 1, 

(3.2) lim W x {f, 0(X O )) = 0, cos(/3)/(l) - sin(/3)/'(l) = 0. 

x— fO 

Here Ao G K and (3 G [0, tt) are fixed, and we can assume that 0(z, x) satisfies 
W(0(\q),9(z)) = (cf. [29, Appendix A]). In this case 9(z,x) continues to satisfy 
(l2~k 

The singular Weyl function is then defined such that 

(3.3) iji{z,x) = (f>(z,x) — rh(z)9(z,x) 

satisfies the boundary condition at x = 1. In this case to becomes the classical m- 
function (see, e.g., [29, Appendix A]) and hence it will be a Herglotz-Nevanlinna 
function satisfying 

1 A 

^A~ ~ TTa 2 

where the spectral measure p satisfies J K dp{\) = oo and J R j^rp < oo. 

Let r G L 1 (0,a) be a positive function on (0,a) and consider the Krein string 
operator 

1 d 2 

r(£) d£ 2 

acting in the weighted Hilbert space L 2 (0, a) and subject to the boundary conditions 

(3.6) 2/(0) = cos0)y(a) - sin(/3)y'(a) = 0. 
Let also c(z,£) and s(z, £) be the fundamental solutions of 

(3.7) -y" = zr(0y, f6[0,o], 



(3.4) m(z) = Re(m(i)) + / ( T — - -— 5 ) dp(X), 
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such that 

(3.8) c(z,Q) = s'(z ) 0) = l, c'0,0) = s(z,0) = 0. 
Define the corresponding m-function M : C \ R — > C as follows: 

(3.9) Vp(z,0 ■= s(z,£) - M(z)c(z,£) satisfies dSU) at C = a. 

The main aim of this section is to establish a connection between the operators 
L and H and hence between the m-coefficients rh and M. 

Without loss of generality we can assume that H > el with some e > 0. Hence 
the solution 8 (x) := 6(0, x) is positive on (0, 1]. Then (see, e.g., [Ml §14]), we set 



(3.10) 



'0 

and define the map U : L 2 (0, 1) — » L^(0, a) as follows 
(3.11) U : w(a;) -> u(£ := ^^(a:). 

First of all, notice that r € i 1 (0, a). Indeed, using (|2.3f) . we get 

r(0«= fto^r= [ 6 2 (x)dx <O c. 
o Jo "oW Jo 

Hence the above definition is correct. Moreover, t7 is isometric: 

ll^ll| 2 = [ a \u((;)\ 2 r(Zn= f TO 2 ^)^T= /V^MIi- 
Jo Jo y oW Jo 

Furthermore, it is not difficult to check that y := Uy solves (|3.7p if y is a solution 
of (12. 2[) . Indeed, this is immediate from the following representation of (|2.2[) (cf. 

Finally, let us show that c(z,£) = U6(z) and s(z, £) = U<j)(z). Notice that 

, r , w-v ,. tffos) /(2Z + l)lim 4e -rt)ffi , y(« > x), Ze(-±,±), 
(£fy)(0) = lim = <^ ( x) / 



and 



1^ 



lim 6>g (a) ^-^-fr = Um „ - ^o) = lim W x (y,6 ) 



5=o x->o dx 6 (x) K-J-0 v 7 x^O 

Hence, using the representation of and 6> from 3 §3], we easily compute that 

(lty)(0) = 0, (C/^(0) = 1, 

and 

(tr«)(o) = i, (uey 6 (o) = o. 

Since both [/(/> and U6 solve (|3.7I) . we are done. Similarly one computes 
(3.12) cos0)(Utl))(a) -sia0)(Utl)Y(a) = 0, cot0) = O (1) ( cot(/3) + l). 
Thus we proved the following result. 
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Lemma 3.1. Let I G [—1/2,1/2) and suppose (|1.2j) . Let also a, r(-), and the 
unitary map U be defined by (|3.10|) - (13.11[) . Lf /3 and /3 are connected by (|3.12l) . 
then the operators H and L are unitarily equivalent: H — U~ 1 LU . Moreover, the 
fundamental solutions (|2.3j) and (|3.8j) are connected by 

c(z,0 = {U<Kz))(£), s(z,0 = {ue(z))(0, 

and the corresponding m-functions defined by p.3[) and ()3.9[) satisfy 

(3.13) m(z) = M(z), z(£m. 

Corollary 3.2. Let I G [—1/2, 1/2) and let Hr 0c \ be a self-adjoint operator associ- 
ated with fll.lj) on (0, c) and subject to some separated boundary conditions at the 
endpoints and c. Then the eigenvalues A„(c) of H/q <c \ satisfy 

71 C 

(3.14) lim 



Proof. Due to the interlacing properties of eigenvalues, it suffices to prove the claim 
in the case of boundary conditions (|3.2j) . By Lemma T3.ll the spectra of operators 
H(o,c) and L(o,£(c)) coincide. However, by [2U §11.8], the eigenvalues A n (£) (= 
A„(c)) of L satisfy 

n 1 f^ C) ATT, 1 T /^7T dt c 

lim - = - / JrHSdt = - / Jei ft) 



□ 

Remark 3.3. ^ Using commutation methods (see, e.g., [30 (|3.14l) can 6e es- 

tended to the case I > —1/2. Let us also mention that detailed eigenvalue asymp- 
totics were obtained in [28] (see also pQj. 

(ii) Note that (|3.14j) implies that the solutions (j)(z,x) and 9{z,x) are entire 
functions in z of order 1/2 and of finite type. Moreover, there is a standard iteration 
procedure to construct the solution c(z,x) (see, e.g., [24, §2]/ and then, by using 
the Liouville transform, we obtain the solution 6{z,x). Let us mention that it is 
difficult to construct 9 (z, x) directly by iteration. Indeed, for the iteration scheme to 
converge it is required that q must satisfy an additional assumption. For instance, 
q G i 1 (0, e) is required if I = . 

3.2. Spectral asymptotics in the case / G [—1/2,1/2). Using the connection 
between the m-functions established in Lemma 13.11 we are able to prove the fol- 
lowing result. 

Theorem 3.4. Let I G [—1/2,1/2), suppose q satisfies (| 1 . 2[) . and let H be given 
by (jl.ll) together with the boundary condition (|3.2[) . Then the corresponding m- 
function (|3.3|) satisfies 



(3.15) m(z) = 3— (1 + o(1)) as \z\ -> oo, 

mi(z) 

where mi is given by (I2.11[) . The latter holds uniformly in any nonreal sector in 
C+. 

Proof, (i) The case I G (-1/2, 1/2). By (JOJ) and Lemma A.2 from [31], we get 
£{x) = (21 + l)x 2l+1 Z(x), | G W^[0, 1], |(0) = 1. 
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Hence, £ as a function of x has limit order 21 + 1 and by Lemma |A.3| the inverse 
x = x(£) of ^(x) has the following asymptotic behavior 

^)=(^ T ) 5 ^(1 + «(!)) as C^O. 
Next, using [3TJ Lemma A. 2], we obtain 

fi r x T l-2l pa 

9 o(t)dt = ————-(I + (i)) = i-(l + o(l)), 



tf(0 = [ r(Oc%= f el 

Jo Jo 



(21 + 1) 2 (1- 21) v w/ A c 



where 

a = Trl' ^-=(1-20(1+20^- 

Therefore, by Theorem IA.21 the m-function (|3.9[) corresponding to the operator L 
satisfies 

M(z)=KxA^ s {-z)-Th(l + (l)) as |«| -> oo; K v = ^ "l^ V 
1 + a (1 — i/) 1 (1 — z/j 

The latter holds uniformly in any nonreal sector in C+. Further, noting that 



(3.17) £ = G(x):= I tt - 77/777^7777 ' R(0 = P(x):= I t log 2 (t)9 2 (t)dt 



1 + a 2 
after straightforward calculations we get 

,~ 1fi x ^ A i+h ,21 M . na r(l/2 + Q sin(7r(/ + 1)) 
(3-16) ^ +| A Q =2 (2/ + 1) r(1/2 _ = 

and Lemma 13-11 completes the proof. 

(ii) The case I = By (|X5| . 6> (x) = V£log(a;)0o(x), where O € W M [0,1] 
and O (0) = 1. Hence by (|3~TT)| we get 

dt 

tlog 2 Jt)6l(t) 
Set 

t" x rlt — 1 r 2 1 

Go ^ := / 7T^77T = wT' Po ^ := / ilog 2 (^ = T (log 2 (x)-log(x)+-). 
Jo tlog^(i) log(or) Jo 2 2 

Note that P, Pq, G, and Go are absolutely continuous and strictly increasing. Let 
P, po, g, and 50 be the corresponding inverses. 

Firstly, since 0q & W 1 ' 1 and 9o(0) — 1, we observe that 

f 9 °( x) dt f x dt 

(3.18) {Go 9o )(x)= / 5^ = / x, x^O. 

y0A Jo tlog 2 (t)e 2 (t) Jo (9 2 og )(t) 

Furthermore, by L'Hopital's rule, P ~ Po and G ~ Go as a; — > 0. Moreover, 

Po{x + o(x)) a log 2 (x) - log(z)(l + o(l)) + 1/2 + o(l) 

— = (l + o(ll) ^ ~ 1, x — ► 0, 

P)(z) V W; log 2 (a:)-log(x) + l/2 

and hence 

(3.19) P{x + o(x))~ P( x ), x^O. 

Noting that R(£) = (Pog)(£) and R (0 = (P °5o)(£) and using (gHH]) and ([3T9]> , 
we get 

(3.20) hm = lim £ ° ^ = lim i^^M = 1. 
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The function Rq, given by 

(3.21) R (0 = (P o g )(0 = 2 + ^ ~ Bo® := ^ 0, 

has limit order oo at and hence, by (13.201) , so has R. By Theorem lA.2[ the asymp- 
totic behavior of m(z) as z — > oo is determined by the function / = F , where 
F(£) :— ^r^y- By Lemma |A.3[ the inverse function / = F~ x is asymptotically 

equal to the inverse / of 

F (0 := — = 2ee 2/ «. 

Note that the function Fo(£) is a bijection from (0,2) to (4e, oo) whose inverse is 
given by 

where W OT is the second branch of the Lambert W- function [35l §4.13] which satisfies 

(3.22) - W m (--) = log(x) + log(log(x)) + 0( l0 f lQg( f )} ), x^+oo 

x iog(x) ' 

Finally, applying Theorem IA.21 we get 

m(pe^) = / (p)(l + o(l)) = -4^(1 + o(l)), P -> oo. 
Lemma |3~T1 completes the proof. □ 

Observe that the m-functions (|2.4|) and (|3.3|) are connected by m(z) = — l/rh(z). 
Therefore, the asymptotic formula (|3.15j) proves Theorem 12.11 in the case I G 
[-1/2,1/2). 

3.3. Spectral asymptotics in the case I > \. Finally the general case can be 
reduced to the previous one as in the proof of Corollary 3.11 from [30]. In fact, 
given I = Iq + k with lo G [— |, h), k G N, we can perform a single commutation 
step to obtain an operator H associated with I — 1 and a new potential q in the 
same class. Now taking into account that <f>(z,x) = (21 + l)x (1 + o(l)) does not 
satisfy our normalization (I2.3P [3H1 Thm. 3.7] implies m(z) = (21 + l) 2 (z — X)rh(z) 
and the claim follows by induction on k. 

Remark 3.5. As pointed in our introduction the solution 9(z,x) is only unique up 
to a term f(z)4>(z,x). Following |31j we will call 9(z,x) a Frobenius solution if 

(3.23) i^-(m)___^ x)s0j 

where n\ := [I + 1/2J . This will fix f(z) up to a polynomial of degree m and ensure 
that m(z) is in the generalized Nevalinna class , where m = [5 + | J (we refer 
to |31] for further details). 

It is easy to see that 9 in Section \3.1\ is a Frobenius solution. Moreover, the 
property of a singular solution 9 to be Frobenius is invariant under a single com- 
mutation. Therefore, it follows from Remark \3.3\ that a Frobenius solution 9 is an 
entire function in z of growth order 1/2 and of finite type. 
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4. Uniqueness results for the inverse spectral problem 

As it was discussed in Remark l3.5l if we choose 6{z, x) to be a Frobenius solution, 
then the m(z) is in the generalized Nevalinna class N£°. This fixes m(z) up to a 
polynomial of degree ni and thus I can be read off from the asymptotics of m(z) in 
this case. 

This observation enables us to show the following improvement of Theorem 8.5 
from [29]. 

Theorem 4.1. Let H(lj,qj), j 6 {1,2}, be the operators defined in L 2 (Q,bj) by 
(jl.ip . Assume that qj satisfies (|1.2[) . lj € [— 5,00), and < bj < +00. Suppose mj 
are defined via solutions 9j(z,x), <f>j{z,x) which are entire functions in z of order 
less than 1. 

If for some c G (0, min{6i, 62}) there is a real entire function g such that for 
every e > 

(4.1) mi(z) - m 2 (z) = g{z) + ( c -2(c-e)|im(y?>|) 

as z — > 00 along some nonreal ray, then l\ = I2 and qi(x) = q 2 {%) for a. a. a; € (0, c). 

Proof. Without loss of generality we can assume that 9j are Frobenius solutions 
such that the m-functions rrij are generalized Nevanlinna functions belonging to 
the class Kj = L2" + fJ- Moreover, (|4.ip implies that g is bounded by a 

polynomial of degree at most n := max{Ki, k 2 } + 1 on two rays in C and hence, by 
the Phragmcn-Lindelof theorem, g is a polynomial of degree at most n. 

Therefore, without loss of generality we can absorb g(z) in the 6j with the larger 
lj and assume g{z) = 0. But then we can read off lj from the asymptotic behavior 
of mj(z) implying l\ = l 2 - But this shows 4>i(z,x) and 4>2{z,x) have the same 
asymptotic behavior and the result follows from [29j Thm. 8.5]. □ 

Corollary 4.2. Let Hj, j 6 {1,2}, satisfy the same assumptions as in Theorem 
\4-l\ Let also the singular m-functions mj be defined via Frobenius solutions. 

If for some c € (0, min{&i, 62}) there is a real entire function g such that (|4.1|) 
holds for every e > as z —¥ 00 along some nonreal ray, then l\ = I2 and qi(x) = 
Q2(x) for a. a. x g (0,c). 

Proof. It suffices to notice that according to Remark |3 . 51 Frobenius solutions are of 
order 1/2 and then to apply Theorem 14. II □ 

Since m{z) is determined by its spectral measure p{\) up to an entire function 
( [29] Thm. 4.1]) we also obtain: 

Corollary 4.3. Let Hj, j 6 {1, 2}, satisfy the same assumptions as in Theorem \4-l\ 
Let also m\ and m 2 be some singular m-functions and p\, p 2 be the corresponding 
spectral functions. 

If there is a real entire function g such that 

(4.2) m x {z) = m 2 (z) + g(z), 

or equivalently p\ = p 2 , then lx = l 2 , b\ = b 2 and qi(x) = q 2 {x) for a.e. x € (0, b\). 

Proof. Without loss of generality we can assume that mj {z) are defined via Frobe- 
nius solutions (this will only change g(z)). But then mj(z) are generalized Nevan- 
linna functions and g(z) can be at most a polynomial. Moreover, we can even absorb 
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g(z) in the 9j with the larger lj and assume g(z) — without loss of generality. 
Hence the result follows from the previous corollary. □ 

Using a different approach based on the theory of de Branges spaces, Corollary 
14.31 was established by Eckhardt |S]. 

As another consequence we also obtain a generalization of item (ii) of Theorem 
2.8 from [25]. 

Corollary 4.4. Suppose H has purely discrete spectrum o~{H) = {A n }. Then the 
eigenvalues X n together with the norming constants 

(4.3) xT 1 = / H*n,x) 2 dx 

Jo 

uniquely determine q and I. 

Proof. This is immediate since in this case the spectrum is purely discrete and 
the spectral measure is uniquely determined by its jumps together with the jump 
heights p({X n }) = In- □ 

5. Quantum scattering theory 

In this section we want to look at the case where in addition to (II .2[) the potential 
has the form 

(5.1) q(x) = - + q(x), q(x) G L\l, oo), 7 G M. 

x 

5.1. The Jost function. Recall that by Weidmann's theorem ([371 Thm. 9.38]) the 
spectrum of H is purely absolutely continuous on (0, 00) with an at most countable 
number of eigenvalues A„ G (— 00, 0], 

Then it is easy to show (cf. [55]) that for k 7^ there exists a unique so called 
Jost solution of rf = zf satisfying the asymptotic normalization 

(5.2) f(k,x) = e ikx -^ los{x \l + o(l)) 

as x — > 00. Here k = *J~z with the branch cut along the positive real axis such that 
< arg(fc) < 2-7T. The Jost solution is analytic in the upper half plane and can be 
continuously extended to the real axis away from k = 0. We can extend it to the 
lower half plane by setting f(k, x) — f(—k, x) = f(k*, x)* for Im(fc) < 0. 
Its derivative satisfies 

(5.3) f'(k,x) = e ite -a M«) (1 + (1)). 

For k G 1R\ {0} we obtain two solutions f(k, x) and f(—k, x) = f(k, x)* of the same 
equation whose Wronskian is given by 

(5.4) W(f(-k,.),f(k,.)) = 2ik. 
The Jost function is defined as 

(5.5) f(k) = W(f(k,.)^(k 2 ,.)) 
and we also set 

(5.6) g(k) = W(f(k,.),9(k 2 ,.)) 
such that 

(5.7) f(k, x) = f(k)6(k 2 ,x) - g(k)cb(k 2 ,x) = f(k)^(k 2 ,x). 
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In particular, the Weyl m-function (|2.4[) is given by 



(5.8) m{ e ) = -9_i^ keC 

/(«) 

Note that both f(k) and g{k) are analytic in the upper half plane and f(k) has 
simple zeros at in n = \/~\^ £ C+ . 

Since f(k,x)* = f(-k,x) for k e K \ {0}, we obtain f(k)* = f(-k), g(k)* = 
g(—k). Moreover, (|5.4I) shows 

(5.9) ( f > ( k 2 ,x) = l^f(k,x)-^f(-k,x), keR\{0}, 



and by (|5 . T[) we get 

(5.10) 2ilm(/(%(fc)*) = /(%(fc)* - f(k)*g(k) = W(f(-k, .), /(fc, .)) = 2ifc. 
Hence 

(5.11) I mM ^)) = - ^g^M = ^ F , fceK \ { 0}, 
implying 

(5.12) dp(A) = X(o.oo)(A) ■ TT^r dA + Y>ndfl(A - An), 



l 

2, 



where 

(5.13) 4>{K,-i 

are the usual norming constants. Since — 7„ equals the residue of m(z) at X n we 
obtain 

(5.14) /(iK») = -2in n 9 ^ Kn \ f(in n ,x) = g(iK„)<j>(\ n ,x). 

In 

By Theorem ED and ([5TT2"1) . 

(5.15) |/(A;)| = 1^(1 + 0(1)), fc->oo. 

As an immediate consequence of Corollary 14.31 we obtain 

Theorem 5.1. Suppose the potential q is locally integrable and satisfies (|1.2p near 
and (|5.ip near oo . TTien i/ie absolute value of the Jost function (|5.5|) on i/ie reaZ 
line together with the norming constants (|5. 13[) and the eigenvalues determines the 
value of the angular momentum I and the potential q uniquely. 

5.2. The phase shift. In this subsection, we assume that the potential q belongs 
to the Marchenko class, i.e., in addition to (|1.2p . q also satisfies 



(5.16) / x\q(x)\dx < oo, q(x) = 



q(x), l >-\, 
\og{x)q(x), l = -\- 

The Jost solution of r/ = k 2 f is defined to satisfy the following asymptotic nor- 
malization 

(5.17) /(fc,x)=e i ^-^)(l + (l)) 
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as x — > oo. Again k — ^fz with the branch cut along the positive real axis such that 
< arg(fc) < 2ir. By the previous considerations, the Jost solution is analytic in the 
upper half plane and can be continuously extended to the real axis away from k = 0. 
We can extend it to the lower half plane by setting f(k,x) — f(~k,x) — f(k*,x)* 
for Im(fc) < 0. 

The Jost function /(.) is defined by ((575]). Notice that f(-k) = f{k)* for k e 
R\{0}. Set f(k) = \f(k)\e-' lS< - k \ k e R\{0}. Note that S(-k) = -S(k) and hence, 
by (15.9[) . we obtain 



(5.18) <f>(k A ,x) = i± ^ Zi siii I kx- — +5(k) I (l + o(l)), x -)• oo. 

The function <5 : R — )• R is called the phase shift and the scattering matrix is defined 

by 

(5.19) S(k) = 4ttt^ = e 2i5 ^, fceR\{0}. 

/O) 

Remark 5.2. In i/ie case g(a;) = we /iaue 

(5.20) ^ x) = d^J?fj l+h ^ x l 

(5.21) f l (k,x)=i x f^HH\ /2 (kx). 



Therefore, 



k 



-I 



(5.22) /,(*,) = lmfc>0 
and hence 

(5.23) 5i(fc) = l, <5/(fc) = 0, fceR\{0}. 

Next we are going to show that f(k) is uniquely determined by its phase 8{k). 
Firstly, observe that the Jost solution (|5.17l) has the form 

(5.24) f(k, x) = C^k- l ^(k, x), Im k > 0, 
where "0 is given by (|B.17[) . Therefore, by Lemma lB.51 the function 

(5.25) F(k) = Cik l f(k) 

is an analytic function in Imfc > satisfying (IB.27[) . By (IB.26[) it extends contin- 
uously to the boundary if / > — i and the same is true in the critical case I = — \ 
except for a possible logarithmic singularity at 0. 
Moreover, (|B.25|) implies 

(5.26) / f^l. dk < oo. 



-1 + 1*1 

Indeed, d(.) is bounded on R and, by (|B.27I) . near infinity it behaves like ImF. 
Hence combining 

C 2 

Imi;i(k 2 ,x) = c/)i(k 2 ,x)lmmi(k 2 ) = -!-k 2l+1 <f)i{k 2 , x). 
with (|B.4|) . (|B.25|) and Fubini's theorem, we obtain the claimed integrability 

ImF(k) M /' 00 ^ i „^m r {kx) 2l x 

(1 + kx) 



dk<C I x\q(x)\ I n \ ,'^ 2l+2 dkdx < C \ x\q{x)\dx. 
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Since q satisfies (|1.2j) and (|5.16|) . the Bargmann bound (see, e.g., [4]) implies that 
H has only a finite number of negative eigenvalues N. Therefore, 



N , 9 \ -1 



is in the Hardy space H°°(C+) and (cf. [27]) we arrive at the following representation 
of the Jost function 

Now Theorem 15.11 immediately implies the following result. 



Theorem 5.3. Suppose the potential q is locally integrable and satisfies (|1.2p near 
cmrf (|5.16p near oo. TTien i/ie phase shift (|5.18[) together with the value of the 
angular momentum I, the norming constants, and the eigenvalues determine the 
potential q uniquely. 

Appendix A. Asymptotic estimates for m-FUNCTioNS 

In this appendix we collect some required results from (3] . Consider the Sturm- 
Liouville differential expression 

-y" + q{x)y = zr(x)y, x 6(0,6), 

where q, r £ Lj oc (0, 6) are real valued and r > a.e. on (0, b). We assume that the 
endpoint x = is regular, i.e., q, r £ L 1 (0, c) for any c £ (0,6). Let m(-) be the 
m-functions corresponding to the Neumann boundary condition at x — such that 

(A.l) s(z, x) + m(z)c(z, x) 

is square integrable near 6 and satisfies a given boundary condition at 6 in the limit 
circle case. 

Define the following function 



(A.2) R(x) := / r(t)dt 

Jo 

and note that R £ W 1 ' 1 ^, c) is positive and strictly increasing on (0, 6). 

Definition A.l. We will say that R has limit order a £ [0, oo) at x — if for 

positive s 

R has a limit order oo at x = if the limit in (|A.3j) equals oo for all s > 1 . 
Define also the function / as the inverse of 

F(x) ' 



xR(x) 

Note that f(y) — > as y — > +oo. The following result is a particular case of 
Theorem 4.1 from [3] (see also [13 Theorem 2]). 
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Theorem A. 2 ([3, 25 ). Assume that R has limit order a at x = 0. Then any 
m-function corresponding to the Neumann condition at x — satisfies 

K 

(A.4) m(pp) = r^p/(p)(l + °(l)) as P^oo, 

where 

(A.5) ^ r( "> 



1 + a' (1 - i/)»T(l - v) ' 

TTie estimates holds uniformly for p in any compact set o/C + . Also. K v = 1 for 

!<). i| . 

Since it will in general not be possible to compute / explicitly, we will rely on 
the following known fact which follows upon combining Theorem 1 and Theorem 3 
from [7j: 

Lemma A. 3 {\J\). Let Rq, R be two functions as in (|A.2|) and define correspond- 
ing functions fo, f as above. Suppose Rq has a limit order a G (0, +oo] and 
Um ^0 igg) = I- Then lim^+oo = 1. 

Remark A.4. Theorem \A.8i was first established by I.S. Kac |23] and Y. Kasahara 
[25] for the case a £ (0, oo). The current form of Theorem \A.2\ was found by C. 
Bennewitz [3] . Moreover, the converse statement has been established independently 
in [3] and [25] . 

Appendix B. Some estimates for the spherical Schrodinger equation 

In this appendix we want to describe some properties of the solutions of the 
spherical Schrodinger equation which are crucial for Section [5] 
The first two lemmas contain estimates for the Green function 

(B.l) Gi(z,x,y) = <f>i(z,x)6i(z,y) - <f>i(z,y)di(z,x) 

and the regular solution <fi(z,x) (see, e.g., [28j Lemmas 2.2, A.l, A. 2]). Here 

(B.2) <h(z,x) = C^z-^^J^^x), 

and 

(B.3) _ 

Y l+ i(^x)- ilog(z)J ;+ i(Vix), Z + ieN , 

where Ji+1/2 and ^+1/2 are the usual Bessel and Neumann functions (see [351 
Chapter 10]). We will abbreviate k = y/z, < arg(z) < 2tt. 

Lemma B.l ([28]). For I > — | the following estimates hold: 
(B.4) ^^(^^e^ 

1-2 „ ,/L| < n( E f 1+ l%V JlmfcKs-y) 



9 ; (z,a;) = -C;z 4 



(B.5) \Gt{k\x, y)\<C[ Y —— : ) ( — ^ ) e l imfc l^, y < 
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and 
(B.6) 



d 

—Gi(k 2 ,x,y) 



< C 



' \ + |% ^ Jl mk \(x-y) 

y 



y < x. 



1+ \k\x / 

For the case I ~ — ^ formula (|B.4|) remains valid and one has to replace ()B.5|) and 
(lB~6l) 6y 

(B.7) 

a; y 



(B.8) 



<+r 



G i (k 2 ,x,y) 



< C 



1%' 



< X. 



Lemma B.2 ([IB])- Assume (|1.2[) and sei = i/ ^ > —1/2 and q(x) 
(1 — log(j^))q(.T) i/ / = —1/2. TTien (f>(z,x) satisfies the integral equation 



(B.9) (f)(z,x) = 4>i(z,x)+ Gi{z,x,y)q{y)(j>(z,y)dy. 

Jo 

Moreover, <f> is entire in z for every x > and satisfies the estimate 



(B.10) \4>(k 2 ,x)~^{k 2 ,x)\<C 
The derivative is given by 
(B.ll) 

and satisfies the estimate 

(B.12) \<j ) '(k 2 ,x)-cj ) ' l (k 2 : x)\<C 



1 + |jfe|a: 



Jim fe|as 



yq{y) 
o 1 + 1% 



4>'{z,x) = 4>'[(z,x) + I —Gi{z 1 x 1 y)q{y)cj){z,y)dy 



l + |fc|a 



Jim ft|as 



yq{y) 
o i + l% 



dy. 



Next we need some estimates for the Weyl solution ip(z, x). We begin with some 
basic properties of the unperturbed Bessel equation. 

Lemma B.3. Let ipi(z,x) be the Weyl solution of the unperturbed Bessel equation, 
ipi(z,x) = $i(z,x) +mi(z,x)<f>t(z,x) = iC i (i\/ z i) /+ 5 J —H^ (iy/^zx) 



(B.13) = [Cl k l +^ ] J^Hi 1 + \(kx) 



where mi is given by (|2.11l) and H^i is the Hankel function of the first kind. 

If I > —1/2, then ipi(k 2 ,x) is analytic in Im/c > 0, continuous in lmk > and 

-i 



(B.14) 



mk 2 , x )\<c 



1 + \k\x 



- |Im k\x 



If I = —1/2, then ip_i(k 2 , x) is analytic in Im /c > ? continuous in C+ \ {0} and 



(B.15) 



hA-i(fc 2 ,z)l <C 



1+ \k\x 



1 - log 



1 + \k\x 



- |Im ft | as 
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Moreover, 

(B.16) Vkip_i(k 2 ,x) Vkxlog(kx), kx -> 0. 

Using the standard technique (see, e.g., [H Chapter 1.5]), one can prove the 
following 

Lemma B.4. Assume (|5.16j) . Then there is a solution ip(z,x) satisfying the inte- 
gral equation 

/>oo 

(B.17) ip(z,x)=ipi(z,x)- G i (z,x,y)q(y)'ip(z,y)dy. 

J X 

Moreover, if I > —1/2, then ip(k ,x) is analytic in Imfc > 0, continuous in 
Irak > and it satisfies the estimate 

(B.18) m\x) -Mk\x)\<c( —^rr-) 1 e"l Imfc l * [°° -^^dy. 



l + \k\xj Jx 1 + 1% 

If l = —1/2, then ip(k 2 ,x) is analytic in Imfc > 0, continuous in C + \ {0} and 
(B.19) 

Moreover, 

(B.20) \Vk$(k 2 ,x)\ = 0(- y/\h\xlog(\k\x)), \k\x->0. 

The derivative is given by 

(B.21) ip (z,x) =<ij> l (z,x) - J —Gi{z,x,y)q{y)ip(z,y)dy 

and satisfies the estimate 

(B.22) I^V)-*^)!^^)"'^- "•jfi?^* 

if I > -1/2 and 
(B.23) 

li>'(k 2 ,x)-f'_,(k',x)l < C (l - log e- M jf -M^. 

if I = -1/2. 

Finally, consider the following function 

(B.24) F(fc) := W^fe 2 ,.),^ 2 ,.)), Imfc>0, fc^O. 

Lemma B.5. Assume (|1.2p and (|5.16|) . Then the function F admits the following 
integral representation 



(B.25) F(k) = l+ 4n{k 2 ,x)(j){k\x)q{x)dx = 1+ / ip{k 2 ,x)<j)i{k 2 ,x)q{x)dx, 
Jo Jo 

Moreover, 

(B.26) F(fc) = cm ; (fc 2 ) + F(fc 2 ), cel. 

where mi is given by 1|2.11[) and _F is entire, and 
(B.27) F(fc) = l + o(l) 
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as \k\ — > oo in Irak > 0. 

Proof. To prove the integral representations (|B.25j) . we need to replace <f> and i\) in 
()B.24|) by (|B.9|) and (|B.17j) . respectively, use the asymptotic estimates for <j), and 
Gi, and then take the limits x — > +00 and x — > 0. 

To prove the second claim, we observe that by (|B.4|) and (|B.14I) 

(B.28) \Mk 2 ,x)4,(k*,x)\<C l+ X , i>-i 

(B.29) |*.j(^,,^,,)| SCT ^ i| .(l-l,(J^ | _)), ! = -i, 

which immediately implies (|B.27|) . 

Finally, the integral representation (|B.25[) implies 

/>oo />oo 

F(k) = l+ 9 l (k 2 ,x)(f>(k 2 ,x)q(x)dx + mi(k 2 ) 4>i(k 2 , x)4>{k 2 , x)q(x)dx 
Jo Jo 

and hence the properties of c and F follow from the corresponding properties of 
4>(k,x) and 9i(k 2 ,x), 4>i(k 2 ,x), mi(k 2 ). □ 

Remark B.6. (i) It is immediate from (|B.26|) that F is analytic in Imfc > and 
bounded and continuous in Irak > if I > — 1/2. 

(ii) By choosing q to be a characteristic function in (|B.25j) . one sees that c in 
(|B.26[) is nonzero in general. In particular, if I = —1/2, then 

(B.30) F(k) = c log(-fc 2 ) + F(k 2 ), ceM, 

where F is entire. 
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